#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
const int N=1e5+10;

int d1[N],d2[N],d3[N];
int up[N];
int p1[N];
vector<int> edge[N];
void dfs1(int u,int fa){
	
	d1[u]=d2[u]=d3[u]=0;
	p1[u]=0;
	for(auto& son:edge[u]){
		dfs1(son,u);
		int len=d1[son]+1;
		if(len>d1[u]){
			d3[u]=d2[u];
			d2[u]=d1[u];
			d1[u]=len;
			p1[u]=son;
		}
		else if(len>d2[u]){
			d3[u]=d2[u];
			d2[u]=len;
		}
		else if(len>d3[u]){
			d3[u]=len;
		}
	}
}
void dfs2(int u,int fa){
	
	for(auto& son:edge[u]){
		int maxSonDist=d1[u];
		if(son==p1[u]){
			maxSonDist=d2[u];
		}
		up[son]=1+max(up[u],maxSonDist);
		dfs2(son,u);
	}
}
int main() {
	
	int n;
	cin>>n;
	for(int i=2;i<=n;++i){
		int t;
		cin>>t;
		edge[t].push_back(i);
	}
	dfs1(1,0);
	dfs2(1,0);
	int nRet=0;
	for(int i=1;i<=n;++i){
		vector<int> vecDist({up[i],d1[i],d2[i],d3[i]});
		sort(vecDist.begin(),vecDist.end());
		nRet=max(nRet,2*(vecDist[3]+vecDist[1]+vecDist[2]));
	}
	cout<<nRet<<endl;
	
	return 0;
}
